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Abstract 



Oh ' The Muskat problem models the dynamics of the interface between 

•^^ . two incompressible immiscible fluids with different constant densities. 

In this work we prove three results. First we prove an L 2 (R) maxi- 
mum principle, in the form of a new "log" conservation law (3) which 
is satisfied by the equation (1) for the interface. Our second result is 
a proof of global existence of Lipschitz continuous solutions for initial 
data that satisfy ||/o||l°° < oo and H^a/olli 00 < 1- We take advantage 
of the fact that the bound ||9a;/o||i<» < 1 is propagated by solutions, 
^ ' which grants strong compactness properties in comparison to the log 

conservation law. Lastly, we prove a global existence result for unique 
strong solutions if the initial data is smaller than an explicitly com- 
putable constant, for instance ||/||i < 1/5. Previous results of this sort 
used a small constant Kl which was not explicit [5, 16, 7, 12]. 
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1 Introduction 

The Muskat problem models the dynamics of an interface between two in- 
compressible immiscible fluids with different characteristics, in porous me- 
dia. The phenomena have been described using the experimental Darcy's 
law that is given in two dimensions by the following momentum equation: 

a 

-u = -Vp- g(0,p). 

Here \i is viscosity, k permeability of the isotropic medium, u velocity, p 
pressure, g gravity and p density. Saffman and Taylor [15] related this prob- 
lem with the evolution of an interface in a Hele-Shaw cell since both physical 
scenarios can be modeled analogously (see also [5] and reference therein). 
Recently, the well-posedness has been shown without surface tension in [6] 
(for previous work on the topic see [1], [18] and [7]) using arguments that 
rely upon the boundedness properties of the Hilbert transforms associated 
to C 1 ' 7 curves. Precise estimates are obtained with arguments involving 
conformal mappings, the Hopf maximum principle and Harnack inequali- 
ties. The initial data have to satisfy the Rayleigh- Taylor condition initially, 
otherwise the problem has been shown to be ill-posed [16], [7]. With surface 
tension, the initial value problem becomes more regular, and instabilities do 
not appear [11]. The case of more than one free boundary has been treated 
in [10] and [12]. 

In this paper we consider an interface given by fluids of different constant 
densities p l , with the same viscosity and without surface tension. The step 
function p is represented by 



p(x,t) =- • 



p 1 , xen 1 ^), 

p 2 , xen 2 {t) = R 2 \n 1 {t), 



for rr(i) connected regions. As the density p is transported by the flow 

pt + u ■ Vp = 0, 

the free boundary evolves with the two dimensional velocity u = (111,112)- 
The Biot-Savart law recovers u from the vorticity given by ui = d xl U2—d X2 u\, 
via the integral operator 

u(x,t) = V ± A~ 1 ia(x,t). 

Darcy's law then provides the relation u = —d Xl p where p/ n and g are taken 
equal to 1 for the sake of simplicity. Then the velocity field can be obtained 



in terms of the density as follows: 

u(x,t) = PV I K(x-y)p(y,t)dy-l(0,p(x,t)). 

Here the kernel K is of Calderon-Zygmund type: 

1 / X\X2 x\ — J 



K ( x ) - | | 2 

■k \ \x\* 2\x\ z 

(see [17]). As a consequence of p € L°°(IR 2 x ]R + ) it follows that the velocity 
belongs to BMO. Moreover, as K is an even kernel, it has the property 
that the mean of K (in the principal value sense) are zero on hemispheres 
[3], and this yields a bound of the velocity u(x,t) in terms of C 1,7 norms 
(0 < 7 < 1) of the free boundary [10]. 

In order to have a well-posed problem we need to consider initially an 
interface parameterized as a graph of a function with the denser fluid below: 
p 2 > p 1 as in [7]. The interface is characterized as a graph of the function 
(x, f(x, t)). This characterization is preserved by the system and / satisfies 

P 2 ~P 1 mr f , (d x f(x,t) -d x f(x-a,t))a 



fj x ,t) = - — —PV I da 



a2 + (/(x,t)-/(x-a,t))2' (1) 

/(x,0) = /o(z), xGR. 

The above equation can be linearized around the flat solution to find the 
following nonlocal partial differential equation 

f t (x,t) = -£^-Af(x,t), (2) 

/(a,0) = / (a), ael, 

where the operator A is the square root of the Laplacian. This linearization 
shows the parabolic character of the problem in the stable case (p 2 > p 1 ), 
as well as the ill-posedness in the unstable case (p 2 < p 1 ). 

The nonlinear equation (1) is ill-posed in the unstable situation and 
locally well-posed in H (k > 3) for the stable case [7]. Furthermore the 
stable system gives a maximum principle ||/||l°°(£) < ||/'||z, 00 (0)» see [8]; 
decay rates are obtained for the periodic case as: 

||/||L-(t)<||/o||L-e- CT , 

and also for the case on the real line (flat at infinity) as: 

,(*)< ll/o||L °° . 
U " 1 + Ct 



L2 1 



Numerical solutions performed in [9] further indicate a regularizing effect. 
The decay of the slope and the curvature is stronger than the rate of decay 
of the maximum of the difference between / and its mean value. Thus, 
the irregular regions in the graph are rapidly smoothed and the flat regions 
are smoothly bent. It is shown analytically in [8] that, if the initial data 
satisfy H^z/oIIl 00 < 1> then there is a maximum principle that shows that 
this derivative remains in absolute value smaller than 1. 

The three main results we present in this paper are the following: 

1) In Section 2, we prove that a solution of (1) satisfies 

2vr J J R J R V v x-a J ) 

= ll/o||£ 2 - (3) 

Furthermore, we have the the inequality 

[dx f d aln(l + ( l^hI^l ) 2 )<C\\f\\ Ll ( S ). 

Jr Jr V v x-a J J 

This identity shows a major difference with the linear equation (2) where 
the evolution of the I? norm provides a gain of half derivative for p 2 > p 1 : 

\\f\\h(t) + (p'-P 1 ) f ds WA^fWUs) = H/olli,, (4) 

Jo 

or equivalently 

ii/ie.w+^r* \*\** ( /(x -' ) " /(g -' ) )' = iiAiii,. 

^ h Jr. Jr \ x ~ a ) 

Notice that this linear energy balance (4) directly implies compactness, 
whereas compactness does not follow from the nonlinear I? energy (3). 

2) In Section 4 we prove global in time existence of Lipschitz continuous 
solutions in the stable case. We understand the solution of (1) using its 
weak formulation: 

dt / dx Tj t (x,t)f(x,t) + / dx t/(x,0)/o(x) 

o Jr Jr /r\ 

/ dt l dx rj x (x,t) PV j da arctan I ■ - 

Jo JR 2vr 7 r \ x-a 



This equality holds Vr/ G C£°([0, T) x R). For initial data /o satisfying 
||/o||l°° < co and ||9j;/o||l°° < 1 we prove that there exists a solution of 
(5) that remains in the spaces f(x,t) G C([0,T] x R) nL°°([0,T]; W ll00 (IR)) 
for any T > 0. We point out that, because of the condition / G L°°(R), 
the nonlinear term in (5) has to be understood as a principal value for the 
integral of two functions, one in H 1 and the other in BMO [17]. 

There are several results of global existence for small initial data (small 
compared to 1 or e < 1) in several norms (more regular than Lipschitz) 
[5, 16, 7, 12] taking advantage of the parabolic character of the equation for 
small initial data. Here we show that we just need Hc^/oIIl 00 < 1> therefore 

fo{x) - f (a) 



< 1. 

a 

Notice that considering the first order term in the Taylor series of ln(l + y 2 ) 
for \y\ < 1, then the identity (3) becames (4). 

3) Our third result, discussed in Section 3, proves global existence of 
unique C([0, T]; H 3 (R)) solutions if initially the norm (6) of /o is controlled 
as || /o ||i < Co where 

ll/ol|i= [ <% \z\\U0\- 

The key point here, in comparison to previous work [5, 16, 7, 12], is that the 
constant cq can be easily explicitly computed. We have checked numerically 
that Co is not that small; it is greater than 1/5. 

2 L 2 maximum principle 

In this section we provide a proof of the identity (3). As we are in the stable 
case, we take without loss of generality (p 2 — p 1 )/(2ir) = 1 to simplify the 
exposition. The contour equation (1) can be written as follows: 

f t (x,t) =PV t c^arctan f /(iM) ~ /(x ~ M) W 

Jr ^ a ) 

We multiply by /, integrate over dx, and use integration by parts to observe 
lci imi2m_ [ A „ f . ., /„a_ + „„ ff(x,t)-f{x-a,t)\ 



o -vr 11/111,2 (*) = - dx da / x (x)arctan 

zat Jr Jr v a 

= — dx da / x (xjarctan 

Jr Jr 



x — a 



We use the splitting 



"&«> -U i 1 ^^-) «*» ( zte ^%- 



1 d 

2dt UJUL * y JrJr^ x — a ) V x — a 

r/«(x)(x-a)-(/(x,t)-/(a,t))\ . //(M)-/(<M)> 



arctan ( - — ]dxda 



x — a / \ x — a 

= h+ h- 

We also use the function G defined by 

rx 

G(x) = x arctan x — In y 1 + x 2 = / dy arctan y. 

Jo 

With these, it is easy to observe that 

h = -l /(s-a)a z Gf /0M) - /(a '' ) ) dxda. 
JrJr v x-a J 

The identity below 

limfr-a)^'*)--^ 

|x|->oo V x — a 

allows us to integrate by parts to obtain 

x — a 



x — a / 



This equality gives 

5SI/&W = -//'" ^( < "'"" fMI )'"^ 

2dt iiit V V x-a / 

and integrating in time we get the desired identity. 

The above equality indicates that for large initial data, the system is not 
parabolic at the level of /. We prove below the inequality 

ln A + ( f( x ^-f( a ^ yy xda < 47rv / 2||/|| Ll (t) 

which shows that there is no gain of derivatives for the stable case. If the 
initial data are positive, then ||/||i,i(£) < ||/o||z,i follows from [8], so that the 
dissipation is bounded in terms of the initial data with zero derivatives. 
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For the proof of the inequality, we denote by J the integral 
J d = f f I In (l + ( /(x) - /(x " a) ) 2 ) dxda. 
We now use that the function ln(l + y 2 ) is increasing to observe that 

JrJr \ a z a 2 J 

The inequality ln(l + a 2 + b 2 ) < ln(l + a 2 ) + ln(l + b 2 ) yields 

J< / [ ln(l + mX f )d X da + [ /lnfl+ 2|/(X : a) ' 2 Wa, 



a 2 / ./in ./lu v a 2 

and therefore 

J<2 f [ In (l + H!Zifli!W a = K 

For if it is easy to get 

K = 2 / dx [ daln(l + 2|/( f |2 

so that an easy integration in a provides 

K = 4ttV2[ dx|/(x)|=47rV2||/|Ui. 

./{x:|/(*)|#0} 

This concludes our discussion of the L 2 maximum principle (3) for (1). 

3 A global existence result for data less than 1/5 

In this section we prove global existence of C([0,T]; i7 3 (]R)) small data so- 
lutions. A key point is to consider the norm 

dei£H/(OI, *>i. (6) 

This norm allows us to use Fourier techniques for small initial data that give 
rise to a global existence result for classical solutions. 



Theorem 3.1. Suppose that initially /o G -ff 3 (M) and ||/o||i < Co, where c$ 
is a constant such that 



2^(2n + l) 2+<5 cl n < 1 



n>l 

for < <5 < 1/2. T/ien there is a unique solution f of (1) that satisfies 
f e C([0, T]; F 3 (]R)) /or any T > 0. 

Remark 3.2. VFe compute the limit case 5 = 0, so that 

2^(2n + l) 2 c !n < 1 
n>\ 

for 



i / 14 x 5 2 / 3 arT — f= — T 

0< c < -i/7 +2y/5(9V39-38) « 0.2199617648835399. 

3 V V 9\/39 - 38 

In particular, 

2]T(2n + l) 2 - 1 c 2 "<l, 
n>l 

»/ sa?/ co < 1/5. 

The remainder of this section is devoted to the proof of Theorem 3.1. 
The contour equation for the stable Muskat problem (1) can be written as 

f t (x,t) = -p(Af + T(f)), (7) 



where we recall that p = p 9 P > and we have 



p'-p 1 



T{!) = i I a./w -«./(*-») ( m ^T da , (8) 

We define 

a 
We consider the evolution of the norm ||/||i (6): 

d 



,,|i(*)= / dei^|sgn(/(0)/t(0 



p/deieiaga(/(0)(-iei/(0-^(r)(0)- 



We will show that the first term dominates the second term if initially 



ll/olli < V( 4 ~^ 13 )/ 6 ' wherey(4-Vl3)/6>l/4. 

The key point, again, is that the constant is given explicitly. 

Notice that under the local existence theorem of [7], this bound will be 
propagated for a short time. Then we may use the Taylor expansion 



°° hl x 2n 



71=1 



to obtain 



Til) = — D" 1 )" / d -^f) i^f? n da. 



(9) 



Notice that 



F{A a f) = /(£)"»(£,«), Hd x A a f) = -i£/(0"i(£,a), 

1 _ e -*c« 

m(£,a) = • 

a 

Therefore 

T(d x (A a f) (A a f) 2n ) = {(-iUm) * (fm) *■■■* (/ro))(£,a), 
with In convolutions, one with —i£fm and 2ra — 1 with fm. Using (9) 

•F(T)(0 = ~^(-l) n / da i '<%!■■■ [ dZ 2n (Z-Z 1 )fe-t 1 )m(Z-Z 1 ,a) 

^ n>1 Jr Jr Jr 

= V / dCl ■ ■ ■ / ^2n(e-6)/(e-6) ( IT /(& -6+l) ) /(6n)M n , 

S 7r ■ /r Vf=i / 

where M n = M n (£, £i, . . . , £ 2 n) is given by 

M n = ^(-l) n /m(e-6,a) I J] m(&-&+i,a)j m(&n,a)da. (10) 
Since m(£, a) = i£/ (is e 4a( - s_1 ^ we obtain 

Mj(^6r--,6n) = m n (£,£l,...,6n)(£l ~ 6) ' " ' (6ra-l ~ 6ra)6n, 



with 



1 /■! 



m n = — I ds\ ■ ■ ■ I ds 2n / da 

7T Jo Jo Jr a 

2n-l 
x exp (ia ^ (s.,- - 1)(£,- - £ j+1 ) + ia(s 2n - l)6n 

3=1 

= - dsi--- / d^nf^ / exp(ia^)— - PV \ exp(iaB) — 

n Jo Jo v Jr a Vr « 

= -/ dsi--- ds 2n (sgnA-sgnB), 
Jo Jo 

and 

2n-l 2n 2n-l 

A = ^2 ( S J ~ 1 )(^' ~ ^ +1 ) + ( S2n ~ 1 ^ 2n = ~^ + 5Z S 5'C? - X] s J'C?'+i- 
i=i i=i i=i 

Additionally 

2n-l 

B = -(£ - 6) + Y, (*i ~ X )^ " &+i) + ( S2 " " X )^n 
i=i 

2n 2n-l 

i=i j=i 

It follows that 

^(r)(0 = E I <%1~' [ dZ2nm n (tel,---,&n)(Z-Zl)fc-tl) 

(2n- 1 \ 

II (^ -^+l)/fe -&+l) ) 6ri(6n), 

with |m n (£,6> . . . , 6n)| < 2. We then have 

/ <% \Z\\nT)(0\ < 2 V / da f d4 1 • • • / d6n ICIIC - 6ll/(£ - 6)1 

JR n>l^ R ^ R ^ R 

X |6 - 611/(6 - 6)1 • • • I6n-1 - 6n||/(6n-l - 6n)||6n||/(6n)|. 

The inequality |£| < |£ - 6| + 16 - 61 + ■ ■ ■ + I6n-i - 6n| + |6n| yields 

m\HT)(o\ <2^(2n+i)( /deiei a i/(oi)( [dm\f(o\ 



2n 

2n+i)u dmv(o\){ijm\no\; 

n>l Jr 
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and therefore 

[dm\HT)(o\ < ( /^iei 2 i/(e)i)2^(2n+ i)ii/n?" 



<( f «wi/(of 



2,^,\2||/||j(3-||/||j) 

2\2 



(I 



Notice ^l^jp < I if < x < \J ^f^- « 0.256400964. If ||/ ||i < 
/- g , then this inequality will continue to hold for some time so that 

|ll/lli(*)<0, 



and we conclude that ||/||i(i) < ||/o||i if ||/o||i < 



4-V13 



Now we repeat the argument but with s > l in (6). Our goal is to obtain 



^||/||2+*(t)<0, 0<6<l/2. (II) 



Let us point out that 



fo||2 + «5<C(||/o|| L 2 + ||ct/o|| L2 ; 



for < 5 < 1/2. Using the inequality 

iei 2+5 <(2n+i)^(ie-^i 2+5 +iei-e2i 2+5 +---+ie2n-i-e2n.i 2+5 +ie2ni 2+5 ), 

we proceed as before to get 

\2+<5m f \\2n 



\t\ 2+S \T(T)(0\dt < / |e| 3+5 |/(OI^2^(2n + l) 



n>l 



In particular, taking ||/||i small enough we find 

/ \c\ 2+5 \nT)(o\dc < I leM/coi*. 

Jm. Jm. 

and bound (ll) therefore holds. 

If ||/||c2,* remains bounded (0 < 5 < l), then from previous work [7], 
we can conclude that there is global existence in C([0,T];H 3 (M)) for any 

T > 0. Since for 

i | \g(x + y) -g{x)\ 
\g\ C 8 = sup j-j^ , 

y^o \y\ d 
ll 



we find 

\g(x + y) - g{x)\ 

\y\ s 

and therefore 



C_ 



\!J 



Jm Jm 



c*,s < c (\\fh~ + f <a ien/(e)i + / di \i\ 2+s \f{o\ 

V Jm Jm 

We conclude that the solution can be continued for all time if ||/o||i is 
initially smaller than a computable constant cq and ||/o||2+<5 is bounded. 
The constant cq defined by the condition 

2^(2n + l) 2+<5 cg n <l, 

n>l 

which has been numerically verified to be no smaller than say 1/5. 

4 Global existence for initial data smaller than 1 

We prove now the existence of a weak solution of the system (1) which can 
be written as follows: 

h = P -d x PV f arctan ( 7(*) ~ / (* ~ a) \ ^ (12) 

where p = (p 2 — p 1 )/2. We first extend the sense of the contour equation 
with a weak formulation: for any rj(x, t) G C^°([0, T) x R), a weak solution / 
should satisfy (5). We show here that this is the case if ||&e/o||l°° < 1- Then 
it follows that ||/[|l°°(£) < ||/o||i°° and \\d x f\\ L oo(t) < \\d x f \\ L °o < 1 as in 
[8]. Then the solution is in fact Lipschitz continuous by Morrey's inequality. 
The main result we prove below is the following: 

Theorem 4.1. Suppose that ||/o||l°° < °° an d \\d x fo\\L°° < l.Then there 
exists a global in time weak solution of (5) that satisfies 

f(x,t) eC([0,T] xM)nL oo ([0,T];^ 1 ' oo (E)). 

In particular f is Lipschitz continuous. 

The rest of this section is devoted to the proof of Theorem 4.1. The 
first step is to prove global in time existence of classical solutions to the 
regularized model (13) below. This is done in Section 4.2. Prior to that, 

12 



in Section 4.1 we prove some necessary a priori bounds. Then, in Section 
4.3 we explain how to approximate the initial data. Section 4.4 shows how 
to prove the existence of solutions of (5), subject to the strong convergence 
established in Section 4.5. 

From now, in the next two subsections we write f = f £ for the solution 
to (13) for the sake of simplicity of notation. The regularized model is given 
by 

f t (x,t) = -eCA^f + ef^ + ?-d x PV [ da arctan(A£/(z)). (13) 

where C > is an universal constant fixed below, and we define 

def f(x) - f(x - a) 



Km 



[a 



with 4>(a) = 4> £ (a) = a/\a\ £ and e is small enough. Initially we consider 
the data /o E W 1 '°°(M.) with ||9x/o|U<»(R) < 1- We will explain how to 
approximate this initial data later on in Section 4.3. 

4.1 A priori bounds 

For the regularized system (13) we obtain the following two a priori bounds 

ll/llz°° < II/oIIl°°, II^b/IU°° < II^x/oI|l°° < i- 

In order to prove the first one, we check the evolution of 

M(t) = max f(x,t) = f(x t ,t). 

X 

Then, where M is differentiable 

M'(t) = f t (x t ,t) 

= -eCK l - £ f{x t ) + ef xx {xt) + ~d x PV [ arctan(A^/(x t ))da. 

We have to deal with the third term above, as the first and the second ones 
have the correct sign. Now 

I(x) = d x PV / arctan(A^/(x))da = d x PV / arctan(A^_ Q /(x))(ia, (14) 

Js. it 
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and thus 

Therefore I(x t ) < (since d x f(x t ) = 0). Then M'(t) < for a.e. t G (0,T] 
and M(t) < M(0). Analogously m(t) > m(0). 
Prom (13) and (14) we have 

/it = — GL-A J x + £Jxxx H -fa:- 

7T 

Using (15) we rewrite I(x) = J 1 (x) + J 2 (x) where 

l) ~ 7r" i^-«i 2 - £ i+(A*u/(x)r a ' 

l) " k H 2 " £ l + (A^/(x))^ a ' 

to find 

Jl{x) = f xx (x)PV f - 1 r—r^ 77 ^ 
7 M 0(x - a) 1 + (A|_ Q /(x))^ 

/• /■ ( x ) _ /(*)-/(<*) i 

1 ' A |x-a| 2 - £ 1 + (A|_ a /(x)) 2 

_ py /■ /s(s)(s - a) - (/(x) ~ /(«)) 2A|_ a /(x) 

7r |x-a| 2 - (1 + (A|_ a /(x)) 2 ) 2 

/ g (x)(x-a)-(l- £ )(/(x)-/(a)) 

X I 12— F 

\x — a\ z e 
and we split further J x (x) = K 1 (x) + K 2 (x) + K 3 (x) + -ftT 4 (x) where 

#V) = f xx (x)PV f — *— — - NV , rfa, 

V ; JxxK > J R <f>(a) 1 + (A-/(x))2 

,. r i \ f(x)~f(x—a) -. 

if 2 (x) = ePy / /xW ~ , 9 i — -i__d«, 

v ; a i«i 2 - e i + (a^/(x)) 2 

r f (x) — /( x )-/( x -») 9 

if 3 ( x ) = -pv / - /xV ; , , 9 — s d a 

V ' 7r |a| 2 ~ £ 1 + (A^/(x)) 2 

14 



da, 



and 



# 4 ( x ) = _py /" da 



|a| 2 - (l + (A«/(x))2)2 

x (/ x (x)|a| e - (l-e)A*/(x)). 



For J 2 it is easy to check that 

/zO) ~ fx{x- a) da 



J 2 x (x)=ePV [ 

JR 



ePV 



|o|2-« l + (A^/(x)) 2 

/x(a) ~ / a (x ~ g) 2(A^/(x)) 2 da 
|a|2- (l + (A^/(x)) 2 ) 2 ' 



Next, as we did before, we consider M(t) = max x f x = f x (x t ,t). Where 
M{t) is differentiable it follows 

M'{t) = f xt (x t ,t) = -eCK l ~ £ f x {x t ) + ef xxx {x t ) + -I x (x t ). 

Now we claim that if M(t) < 1 then M'(t) < for a.e.t. We can conclude 
analogously for m(t) = min x f x > —1, m'(t) > for a.e.t. 
We check that if M(t) < 1, then 

-eCA 1 " 6 /^**) + E/aaxCxt) + ^J s (x t ) < 0. 

7T 

We can use the following formulas for the operator A 1 " 6 f x 

^- £ Ux) = Cl (e) / fx(x) ~ ' xi ? ~ a) da, (16) 

Jr |a| 

A 1 " 6 /.^) = c 2 (e) / Jx{ | , 2 £ a da, (17) 

Jr l«l 

where < c m < ci(e), c 2 (e) < cj\/ for < e < 1/4. 
We claim that 

-eCA 1 - 6 /*^) + ^(K 2 (x t ) + J 2 x {x t )) < 0. 

7T 



We will show that 



-e^-k^Uxt) + ^# 2 (x t ) < 0, 

A TT 
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using (17) and that 



by (16). In fact 



Z IT 



-e^A^f x ( Xt )+P K 2 (x t ) 

Z 7T 

7r |a| 2 ~ e 1 + (A^/(x t )) 2 

The mean value theorem gives 

\f{x)-f{x-a)\/\a\<\\f x \\ L ~. 
Thus if we take C > — &- we obtain the first inequality. Also 

-e^A^f x ( Xt )+Pj\ Xt ) = 

Z 7T 

pv f Ux t ) - f x (x t - a) ^(A'/fe))' + g^g ~ f , 
A l«l 2 " £ l + (A*/(z t )) 2 

p „ /" /^Q-Zx^-a) 2(A £ a /(x f )) 2 da 
7r H 2 " e (1 + (A|/(x t )) 2 ) 2 ~ 

We find fxxxixt) < and i^ 1 (x t ) = 0. We still have to deal with K 3 and 
K 4 . Considering K 3 (xt) + K 4 (xt), we realize that if 

P(a) = 2 + 2(AU(x t )) 2 + 2(A%f(x t ))(f x (x t )\a\ £ - (l-e)A* /(x t )) > 0, 

we are done. We rewrite 

P(a) = 2 + 2e(AU(x t )) 2 + 2(A%f(x t ))f x (x t )\a\ £ , 

and therefore we need 

\(AU(x t ))f x (x t )\a\ £ \ < 1. 

This fact holds if 



\f{xt)-f(x t -a)\ 1 



|/| c i-2 e =sup- 
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Now we will check that if ||/||l°° < H/olli 00 and H/zIIl 00 < 1 then |/| c .i_ 2e < 
1 for e small enough uniformly. We replace 2e by e without lost of generality. 
If II/oIIl 00 = or ||/x||l°° = all done. Otherwise 

\f(x t ) - f(x t -a)\ ... , x£ 

— S \\jx\\L°°0 , 



Ia| 1_e 



for < \a\ < 5 and 



\f(x t )-f(x t -a)\ <2 J|/olU° 



gl-e ' 

for \a\ > 5. We take 5 1_£ = 2 1 1 /o 1 1 z,°° / 1 1 /a; 1 1 z,°° and therefore 

|/|oi-« < max{\\f x \\ L <x>, ||/x||^ t ^(2||/ ||l-) t ^}. 

Now it is clear that given ||/o||l°°, if H/xlU 00 < 1 there exists £o > such 
that, for any < e < £o it follows that l/lc 1 - 8 ^ 1- 

4.2 Global existence for the regularized model 

We use here the a priori bounds to show global existence. Local existence 
can be easily proving using the local existence proof for the non-regularized 
Muskat problem (1), as in [7]. We use energy estimates and the Gronwall 
inequality. As we did for (1), it follows that 

|ll/lbW = -^//^ln(l + ( /fet) - /(a ' t> )>^ 

dt tt J R J R \x-a\ e V V x-a J J 

-2C£\\k { - l -^f\\ L2 (t)-2£\\f x \\ L 2(t). 
Therefore ||/|| L2 (t)<||/ || L2 . 
Remark 4.2. The theorem of this section can also be found with 

H/oIIl 2 < °° instead of ||/o||l°° < oo. 

We picked the version above because it is more general. We see that if the 
solution satisfies initially a L 2 bound then f(x,t) G L°°([0,T];L 2 ' 



Next, we consider the evolution of 

/ dlfdlf t dx < -C£\\A^~^ 2 d 3 x f\\ L 2(t) - £\\d x f x \\ 2 L2 + L l + L 2 , 
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where 



L x = P -\ dlf{x)dl(pV f fx{x) j x [ X a) da)dx, 
kJr xK,x ^ Jr </>(«) > 



L2 " -* k dJ{X)dx[PV L W) 1 + (A« /(*))' 

The term f x (x) cancel out in L\ due to the PV and an integration by parts 



shows that 



Li = --C(e) [ dlf(x)K 1 - £ dlf(x)dx < 0. 



For L2 one finds 

Wr V Jr 0(a) 1 + (Ag,/(x))V 

and the splitting L2 = Mi + M2 + M3 gives 

x ~*h xm h *(«) i + (A%f( x )y dx > 



3 P f cAst \ f d xf( x ) ~ d lf{x - a) f x (x) - f x (x - a) 



M 2 Jt f d tf(x) [ 

7T Jr Jr 



(a) 0(a) 

„ 2(A%f(x))da 



;i + (A|/(x))2)2 



dx, 



3 " - L xt[) L V 0(a) J (1 + (A S /(s))*)3 dX - 

For Mi we proceed as follows 



\Mi\ = - I / da dx + da dx) 

n \J\a\>i Jr J\a\<l Jr J 

<C(e)(\\f\\ L ~+l)\\d 3 x f\\ L 4d x f\\ L 2. 
The identity 



Jo 



d x f{x) - d x f(x -a)= d x f(x + (s - l)a)ads, 
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yields 

6p f 1 



\M 2 \<^ f ds f f -^ F f dx\dif{x)\\dlf{x + {s-l)a)\ 

K Jo J\a\<l M £ JR 



+ - t ds I Tl^ i dx l^/(x)||^/(x+(,-l)a)| 

^ Jo J\a\>l H ' iS Jr 



x(\f x (x)\ + \f x (x-a)\) 

6p f 1 , /" da 

/ ds / 

'0 J\a\>l 

x(\f x (x)\ + \f x (x-a)\)(\f(x)\ + \f(x-a)\), 
and therefore 

|M 2 |<C(e)(l + ||/||l«.)||^/||l2||^/|| l2 ||/ x || L cc. 
In M3 we use the splitting M3 = JVi + iVo, where 

N\ = — / do / (ire, N2 = da dx, 

n J\a\>l JR J\a\<l Jr 

and then 

\Ni\ < ^H/.Hi- /'(fa/ ,-^a? /ifa:|^/(x)|(|/«(x)| + |/ x (x-a)|) 

7T Jo J|a|>l l«| J d£ JR 

<<mili~ll/*IM^/ll^ 
^cii/.Hi-di/iiL. + 11^/11^)11^/11^. 

To finish, the equality 

/x(») - fx(x -a) = dlf{x + {s- l)a) a ds, 







allows us to obtain (since + 4 + 4 = ^) : 



2 ' 4 ' 4 

|iV 2 |<— \\fx\\ L oo! dr [ ds f -^3- /dx 
71" Jo Jo J|a|<i H de Jr 

x \dtf(x)\\d 2 J(x + (r - 1)^)11^/(1 + (a - l)a) 

4rll Ila2jrll2 



<c||/ :c ||^ii^/llL 2 ||^/iii4. 

The following estimate 



l^ 2 /H 4 L 4 = [{dlffdlfdx = -3 [ (d 2 J) 2 d x fd x fdx 

Jr Jr 



<3||/,||iH|flg/||i*||flg/|| 



1,2 1 
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yields 

|iV 2 |<C||/ x ||| 00 ||^/|| L2 ||^/|| L2 . 

Using Young's inequality 

J t \\d 3 J\\h < C(e)(||/||t- + ||/||i- + ll/.llt- + H/.Hioo + 1) 

x(ll/ll£> + ll#/lllO. 
and therefore the Gronwall inequality yields 

h(t) + \\d 3 J\\Ut) < (\\fo\\h + \\%fo\\h)e*P (J*C(e)G(s)ds 
for 

G(S) = \\f\\U + ll/llloc + ||/x||!oo( S ) + ||/*||£oo(*) + 1. 

We find / G C([0, T]; fl" 3 (M)) for any T > by the a priori bounds. 

4.3 Approximation of the initial data 

The approximation of the initial data described below is needed in order to 
construct a weak solution. First consider a common approximation to the 
identity ( G C™(R) satisfying 

dx ((x) = 1, C > 0, C(x) = C(-x). 

Now we denote the standard mollifier Ce( x ) = C( x / e )/ £ s0 that Ce( x ) con_ 
tinues to satisfy the normalization condition above. 

For any /q G VF 1,00 (1R) and HSr/ollz, 00 < 1) we define the initial data for 
the regularized system as follows 

/o(x) " 1 + ex* ■ 
Notice that /jf G H S (R) for any s > 0, and 

||/olk°° < ll/olU 00 - 

More importantly, 1 1 ^c /q 1 1 z,°° < H^s/olli 00 if e is sufficiently small (here e 
will generally depend upon the size of ||/o||l°°)- I n particular 
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and clearly 

(Ce*d x f )(x) 



l + ex 2 



< ||(Ce*3r/o)||L°°(R) < ||§c/o| 



On the other hand, by splitting into |x| < e 2 ' 3 and \x\ > e 2 ' 3 we have 
that 

2ex(l + ex 2 )~ 2 < 2max{e 1/3 ,e}. 

On the unbounded region we have 

x{l + ex 2 Y 2 = (— +ex 3/2 J < 1. 

Thus, the desired bound follows if e is small enough. Therefore global exis- 
tence of the regularized system (13) is true with /q if e is small enough. 

Now consider the solution to the regularized system (13) with initial 
data given by the /q just described above. For e > sufficiently small, we 
decompose 



■nix, o)f^(x)dx = / n(x, o) ^ / UA 2 > dx = i{ + if, 



( Cs*fo)(x ) 

1 + ex 2 



where 



and 



q= [ r 1 (x,0)(Ce*fo)(x)(—^-l)dx, 
Jr \ l -\- ex / 



r](x,0)((e* fo)(x)dx. 



We apply the dominated convergence theorem to find that as e \. it holds 



that If — ^ 0. For If we write 



n 



[ ( £ *( V (;0))f (x)dx. 

JR 

The L 1 approximation of the identity property shows that 

If -»• / 77(x,0)/ (x)dx. 
Thus, it remains to check the convergence of the rest of the terms in (5). 
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4.4 Weak solution 

In this section we prove that solutions of the regularized system converge to 
a weak solution satisfying the bounds 

||/IU»(t) < ||/o|U», \\d x f\\ L °°(t) < ||0 X / O |U» < 1. (18) 

Given a collection of regularized solutions {/ e } to (13), we have the uniform 
(in e > 0) bound 

||/ £ ||L-(t)<||/o||L-, ||^/ £ || LOO{R) (t) < 1, e>0. (19) 

This implies that there is a subsequence (denoted again by f £ ) such that 

/ f £ (x,t)g(x,t)dxdt — > / / f(x,t)g(x,t)dxdt, 

o Jr Jo Jr 

/ d x f £ (x,t)g(x,t)dxdt — > / / d x f(x,t)g(x,t)dxdt, 

o Jr Jo Jr 

for / G L oo ([0,T];VF 1 ' oo (IR)) and any g G ^([0,7] x R) by the Banach- 
Alaoglu theorem. We find weak* convergence in L°°([0,T]; W 1,00 (M.)). 

We denote B^ = [—N,N\, then we claim that there is a subsequence 
(denoted again by f £ ) such that 

II/ £ -/Hl°°([o,t]xB jv ) ->0, as e^O. 

We will prove this fact in the next Section 4.5. Then, up to a subsequence, we 
find uniform convergence of f £ to / on compact sets. Since f £ G C([0, T] xM) 
we find that / is continuous. 

The only thing to check is that as e 4 we have 

dt \ dx r)Jx,t)-PV da arctan ( - , , ~ — 

Jr ^ Jr V <P £ ( a ) 

-)• / (it / dx r] x (x,t)-PV da arctan ( - ^ - - ^— - a ' 

Jo Jr ^ Jr \ a 

where 4> £ (a) = a/\a\ £ . The rest of the terms will converge in the usual 
obvious way (since they are linear). 

Choose M > so that supp(?]) C Bm- For any small 5 > and any 
large R^> 1, with R > M + 1 we split the integral as 

da - I da + da + / da. 

Jb r -b 6 
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We begin with a proof that the first and last integrals separately are ar- 
bitrarily small independent of e for R > sufficiently large and for 5 > 
sufficiently small. One finds that 



arctan 



f{x) - f(x - a) 



[a 



< 



TX 



Here we don't need any regularity for f £ , and we conclude 

dt / dx r) x (x,t)—PV / da arctan I r~r~\ 

Jm 7T Jb s V <P\<X) 

< P\\ r lx\\L 1 ([0,T]xR)fi- 

Therefore this term can clearly be chosen arbitrarily small, depending upon 
the smallness of 5. 

We now estimate the term integrated over B^. We note that 



arctan y 
and therefore 



d 



(aict&n(sy))ds = y . 
o d s Jo ± + s 2 y 2 



1 



ds, 



V Jo ! + s y 

This is morally the first order Taylor expansion for arctan with remainder 
in integral form. With this expression we have that 



PV I da arctan 

IB' 



PV 



da 



B', 



r( x ) - r( x - a) 

cj) e (a) 

r( x )-r(x-a) 

(j) E {a) 



-Hur 



s 2 ds 



Here H^ is a (Hilbert-type) Transform, which has the form 

f £ (x - a) 



HUf £ ) = PVl da 

<B% 9 («) 



The principal value is evaluated at infinity (if necessary). For the second 
term in the left hand side notice that the integral is over B^ and the principal 
value is not necessary. In particular, we have 



l 
da I ds 



B>, 







<c\\r\\i 



- CH/oIlL 



R a 



3-3e 



< 



R 
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This term is therefore arbitrarily small if R is chosen sufficiently large. We 
are going to show the same for 



fM 
I R = / dx Vx ( X ,t)H R (f). 
J-M 



We write Ir = Jr + Kr where 

" M f £ (x - a) 



Jr = lim / dx rj x (x,t) I da 



<x(x,t) 

M J-n 9 £ {a) 

/M rn f £ (x — a) 

dx r] x (x,t) / da — — — - — . 
-M Jr <r\pt) 

We shall show how to control Jr, the same follows for Kr. We write 



/M rn 

dx rj x {x,t) / da 
-M Jx+R 



/ £ («) 



4> £ (x — a) 
An integration by parts yields 

,. f M , / n ff{x + R)\R\ £ „ N f n , f £ (a) 

Jr = km / dx V (x, t) ^ -f-L + (1 - e) / da , ,' 

n^ooJ_ M \ -R J x+R \x-a\ 2 £ 

Hence 

|Ji?|<2||r ? || L i||/ ||L-/ J R 1 / 2 . 

Since tke same estimate kolds for \Kr\, one finds tkat Ir is arbitrarily small 
if R is arbitrarily large. 

It remains to prove tke convergence of 

T f , ^ ( f £ {x)- f £ {x-a) 

I da i] x (x,t) arctan r~T^ 

o Jb r -b s V £ (a) 

Recall tkat we kave uniform convergence on compact sets. Lets consider 

f £ (x) - f £ (x-a) 



G e 



cj) £ (a) 



wkere x G Bm and a S Br — B$- Since arctan is a continuous function tken 
arctan(G £ ) — > arctan(G°) converges uniformly. Tkus also tke integral over 
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a bounded region of arctan(G £ ) also converges. Then for any R > M + 1 
and any small 5 > as e \. we have 



9 



dt / dx rj x (x,t)— I da arctan 

Jm. kJbr-Bs 

dt / dx rj x (x,t)— I da arctan 



fix) - F(x - a) 



4> £ (a) 
f(x) - f(x-a) 



o 



We conclude by first choosing R sufficiently large and 5 > sufficiently small 
and then sending e 1 0. Note that R and 5 will generally depend upon the 
size of ||/o||oo) but this has no effect on our argument. 

4.5 Strong convergence in L°°([0,T];L°°(B R )) 

In order to prove the strong convergence in L oo ([0,T];L°°(Bji)), the idea 
is to use the non-standard weak space W* '°°(Bji) which will be defined 
below. Crucially, we will have the uniform bounds: 



sup 

te[o,r] 



sup \\f(t)\\ w i^ (BR) <C||/o|| L o 
te[o,T] 

£ 

< c||/oI|l°°(m), 



df 



(20) 



wr 2 -°°(B R ) 



where C does not depend on R or e. Prom here we will conclude that for 
any finite R > there exists a subsequence such that f e — > f strongly in 
L°°([0,T\;L°°(B R )). 

We define the space W* '°°(Bji) as follows. For v G L co (B R ) we consider 
the norm 



-2,oo 



sup 

4>GW^(B R ):\ 



1<1 



Bt- 



<fi{x)v{x)dx 



Here W ' (Br) is the usual set of functions in W 2,1 (Br) which vanish on 
the boundary of B R together with their first two weak derivatives. Now 
the Banach space W* ,co (Bfi) is defined to be the completion of L°°(B R ) 
with respect to the norm || • ||-2,oo- I R general this is all we need for our 
convergence study. This will be explained after the proof of Lemma 4.3 
below. The full space W*~ '°° may be a large space, but because we are 
going to deal with f £ and ^, both in L°°([0, T};L°°(B R )), it is not difficult 
to find the norms of both functions in W* - '°°(B R ). 
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These spaces are suitable because 



W>°°(B R ) C L°°(B R ) C W-*>°°(B 



R 



-2,oo, 



Now the embedding L°°(Br) C W* z °°(Bji) is continuous, and the embed- 
ding W 1,00 {Br) C L°°{Br) is compact by the Arzela-Ascoli theorem. 

We now proceed to discuss the convergence argument. Arguments re- 
lated to Lemma 4.3 below are described for instance in [4]. However in [4] 
reflexive Banach spaces are used. None of thespaces used here are reflexive. 

Lemma 4.3. Consider a sequence {u m } in C([0,T] x Br) that is uniformly 
bounded in the space L°°([0,T]; W 1,00 (Br)). Assume further that the weak 
derivative -^ is in L co ([0,T];L co (Br)) (not necessarily uniform) and is 



dt 



-2,00, 



uniformly bounded in L°°([0, T]; W* '°°(Br)). Finally suppose that d x u m G 
C([0,T]xBr). Then there exists a subsequence ofu m that converges strongly 
inL™(\Q,T];L™(B R )). 

Proof. Notice that it is enough to prove that the convergence is strong in 
the space L°°([0, T]; W* - '°°(-Br)) because of the following interpolation the- 
orem: for any small rj > there exists C^ > such that 

\\u\\l°° < V\\ u \\l,oo + C T? ||ti||_ 2 ,oo- 

This holds for all u £ W l, °°{Bp;). See, for example, [4, Lemma 8.3]. Here 
we can replace reflexivity with the Banach- Alaoglu theorem in W 1,co (Br). 
Let t,s £ [0, T] be otherwise arbitrary. We have 



u m (t) -u m (s) = 
This holds rigorously in the sense that 

/ u m {t)4>dx - / u m {s)4>dx 
Jb r Jb r 

for any (j> e W 2,1 (Br). Clearly, we find 

\\u m (t) -Um(s)\\ 
and therefore 



1 , du 
dr 



0t 



d.T 



(r). 



du r 



B T - 



0t 



(T)4>dx. (21) 



{ R> re[0,T] 



du r 



Dt 



(r) 



wr 2 '°°(B R ) 



\u m (t) -u m (s)\\ 



W» (B R ) — I 



(22) 
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where 

du 
L = sup sup 

mSN T e[0,T] 



<r) 



W- 2 ^(B R ) 



dr 

Now we consider {t fc } fc6 N = : [0,T] n Q. We have u m (t k ) G M^ 1,00 (-B fl ) for 
any m and fc. By the standard diagonalization argument, we can get a 
subsequence (still denoted by m) such that 

Um(tk) -> «(*fc), 

for any fc in L°°(Br) as in the Arzela-Ascoli theorem. 

Consider e > 0. Since [0, T] is compact, there exists JgN such that 



[0,T]cUk-^+-^ 



6L' J 6L/ 

Then there exists iVj such that V mi, vni > iVj it holds that 

\\u mi (t k:j ) - , u m2 (t fci )|| H/ -2,cx) (BH) < e/3. 
Taking N = max iV,-, V mi, 7712 > iV it is easy to check that 

sup ||u mi (t)-u m2 (t)|| w -2,cx. fB x <e. 
te[o,T] * v R > 

We find the sequence uniformly Cauchy in L°°([0,T]; W* - '°°(-Br)) which 
therefore converges strongly to an element in L°°([0,T]; W*T '°°(-Br)). D 

Now we apply the above Lemma 4.3 to prove the strong convergence 
which was claimed in the previous Section 4.4. It remains only to prove 
that for any solution f £ to (13) we have ^ in L°°([0, T]; L°°(B R )) (but not 
uniformly) and that the second inequality in (20) holds for all sufficiently 
small e > and for any R > 0. 

Recall that f £ G C([0,T];H 3 (R)), then in (13) the linear terms are 
bounded easily. The nonlinear term can be written as 

NL ~ (1 £)CA f L na) l+wwr**' 

and therefore 

\NL(x,t)\<C(e)\\r\\ H3 (t), 

by Sobolev embedding. 
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<9/ e 



2,00/ 



The norm of -^- 6 VF* '°°(Bji) is given by 



<9f e 



= sup 

W- 2 >°°{B R ) <t>&W^\B R ):U\\ w 2,i<l 



I dx -Trr( x i t ) ( t )i y x ) 



Since (ft vanishes on the boundary of Br, we can think of (ft{x) as being zero 
outside of the ball of radius R. Then we are allowed to integrate over the 
whole space M. This is also unimportant below because the functions we 
are estimating are defined on the whole space ]R. It is however important 
because we want to estimate the non-local operator A 1_e in this norm via 
"integration by parts" . Then we have 

I=i A 1 ~ £ f(x)(ft(x)dx = I A 1 - £ f(x)^(x)dx = I f(x)A 1 - £ ^(x)dx, 
Jb r Jr. Jr 



and therefore 



We compute 



III^H/IU-WHA^IUi. 



a 1 -^) = c / ^:t {x r a) d a 



K> 



2-e 



da + da = J\{x) + J2(x), 

\a\>l J\a\<l 



thus 



\J\{x)\dx < 



da 

1 2-e 



'|a|>l \ a 

It is easy to rewrite J2 as follows 



dx(\(ft(x)\ + \(ft(x-a)\)<C\\(ft\\ LHBR) . 



J 2 (x) = c 



(x) — 4>{x — a) — (ft x {x)a 



|a|<l 



n 



2-e 



da, 



and therefore the following identities 



<ft(x) — <ft(x — a) = a I (ft x (x + (s — l)a)ds, 
Jo 

(x) - (ft(x - a) - (ft x {x)a = a 2 / (s — l)ds dr (ft xx (x + r{s - I) a), 



allow us to find 



\Ji(x)\dx < / ds dr I dx\(ft xx {x+r{s-l)a)\<2\\(ft xx \\ L i, BR) . 

J\a\<iJo Jo Jr 
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Now we clearly have using the a priori bounds that 

11^- £ f £ \\w- 2 -°°(B R ) + WfxxWwr 2,00 (B R ) - C H/ £ |li°°(R) - C II/oIIl°°(R) • 

Here the constant is independent of e and R > 0. 

For the last term in (13) we consider the duality relation 

/ dx d x u(x) PV / da arctan(A^,/ e (x)). 

Using exactly the arguments from Section 4.4 with say R = 5 = 1 we have 



dx d x u(x) PV l da arctan(A^/(x)) 



< (711^11^1,1 ||/o||ioo( R ) • 
We thus conclude (20). Q.E.D. 
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